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Abstract 

A new version of application Pauli-Villars regularized Green functions in 
the quantum field theory using higher derivatives is proposed. In this version 
the regularizing mass M is large but finite. Our approach is demonstrated 
and discussed on the example of QED. It is shown that in our case there 
are no ultraviolet divergences and - on the example of the selfenergy spinor 
Feynman diagram - no infrared ones. 
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1 Introduction 



One of the unwritten rules of the standard quantum field theory is the limited 
order of the derivatives in the field equations. As is well known this order is 
two for the field with integer spin and - one for half spin. This rule appears, 
at first, with the Klein-Gordon and Dirac equations, which stand in the origin 
of many field models. We can suppose that the analogy with the classical 
mechanics is one of the main reason for this situation - let us recall that only 
first order time derivatives enter the action of the usual mechanical systems. 
However, mechanical actions with higher time derivatives were considered 
sometimes by different authors ]]J . In this case they use the special canonical 
formalism first developed by Ostrogradski 

Recently models with higher derivatives appeared also in quantum field 
theory, especially in the two dimensional one (see for example [|[ where the 
reader can find other references on this subject). However, we are going 
to concentrate our attention to the use of higher derivatives in the 3+1 di- 
mensional quantum field theory. In this direction we can point out works 
in which the authors study several mathematical properties of the free field 
equations with higher derivatives. These equations are characterized by a 
polynomial in the d'Alambertian operator It is interesting to mention 
the work ||, where the higher derivative kinetic term was introduced in the 
Nambu-Jona-Lasinio model. 

In the last few years appeared several works in which the authors study 
the so called "differential regularization" first proposed by Freedman, John- 
son and Latorre ||. Though the derivatives are used there only for the 
regularization of the loop diagram in the coordinate representation, we may 
suppose that these formal rules maybe follow from the appropriate quantum 
field theory with higher derivatives. 

The main difficulty arising when we use higher derivatives in the quantum 
field theory is the appearance of the indefinite metric in the state space and 
all physical consequences from this. That is why quantum field theory with 
higher derivatives in the free part of the corresponding Lagrangean can not 
exist, if we want every vector from the state space to have a physical meaning 
(here we have in mind the spaces of in- and out- states). However, there exists 
another approach in which the physical states form the appropriate subspace 
with definite metric in the frame of the whole indefinite state space. Such 
an approach was sketched by Hawking 0. In this case the physical theory 
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appears as immersed in a wider theory, analogously to the situation in the 
quantum electrodynamics in the Lorentz gauge. 

Different authors adduce different arguments for the application of higher 
derivatives depending on the concrete problems they pay attention to. The 
presence of higher derivatives in the kinetic part of the field Lagrangean leads 
to free propagation function in which the power of the momentum p is less 
than -2. This fact means that the divergences of the Feynman diagrams 
will become smaller. For example if the propagator of the fermions in the 
usual spinor QED has the behaviour p~ 3 , when p — > oo, then all well known 
divergent diagrams become finite. This fact is the main argument for us 
to consider higher derivatives. The aim of the present paper is to give an 
alternative formulation of the spinor QED in which the fermion field obeys 
a third order differential equation. 



We are going to consider a model in which the main fields are the spinor field 
<p(x) and vector one A^(x). Besides we have an auxiliary spinor and scalar 
fields denoted by ip(x) and <&(#) respectively. To write down the action S of 
our model, let us introduce the following notation 



where e is the dimensionless electric charge. (From here on the Greek indices 
such as ([/,, v, A) run from to 3 - as vector indices - those such as (a, (3, 7) 
run from 1 to 4 - as spinor indices). Then 



where the parts L sp and L v of the Lagrangean have the following form: 



2 The Model 





(i) 






L V (A) 




(3) 
(4) 
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m and M are mass parameters. Here with V* and <p we have denoted the 
complex and Dirac conjugation respectively. Moreover we have denoted: 



uD v = uD^v — D^u.v 



L T ($) is the free kinetic part of our Lagrangean for the scalar field $. As we 
can see below the concrete form of this term does not matter for our model. 

For the construction of the Lagrangean (|3|) we have used, in general, a 
known manipulation for the exclusion of the highest derivative ( in our case 
- third derivative), with the help of the auxiliary field ip{x). It is easy to 
verify, that if we replace the field ip(x) by new a auxiliary field x{ x ) with the 
relation 

ip(x) = x(x) j= = (i^Vf, + m)<p(x) (5) 



m 



then we can obtain the following new form of the spinor part of the La- 
grangean 

2 

+ W- - m 2 ^/27"V, + m)y 

+ ^(V2 7 ^ + m)^ + x(V 2 7^-m)x (6) 

without the change of the action (0). We see that the free field x is n °t 
connected with the rest of the fields and has no significance for our model. 
In spite of the equivalence between the spinor part (^) of the Lagrangean 
and L sp {ip) from (|6|) we give our preferences to the one from eq. @, because 
it is more suitable for the applications of the canonical quantum formalism. 

The Lagrangean parts L ap ((p) and L V (A) are gauge invariant under the 
following local transformations: 

<p(x) — > expierj(x).ip(x); ip(x) — > expierj(x).ip(x) 
A^x) -> A^x) + d^x); -> + rj(x) (7) 

where rj(x) is the corresponding local parameter. 
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Let us try then to write down L T ($>) in the gauge invariant form. In this 
case it has to depend on — A^, because this is the only gauge invariant 
combination in the presence of the field $: 

L T {$) = W(d^ - A,) (8) 

In this case the action (||]) becomes gauge invariant. However, with a simple 
redefinition of the fields ip, tp and A^ with the help of the gauge transforma- 
tions (0), the field $ will disappear from the action. This means that the free 
Lagrangean L T (Q) can not be gauge invariant for arbitrary 7](x). Taking into 
account that the conformal dimension of $(x) is zero, it is natural to suppose 
that the free part of its equation must be D 2 $ (□ is d'Alambertian opera- 
tor). Then the gauge invariance is of the type (|7|) but with r](x) satisfying 
the equation 

U 2 r]{x) = (9) 

As we shall see below this restriction is in accordance with the usually used 
gauge fixing rules in QED. 

Now we can give the explicit form of L T ($): 

L T ($) = d^d„U + hj 2 (10) 

U(x) is an auxiliary field with the help of which L T (<&) depends on first 
derivatives only. The gauge transformation of U (x) has the form 

U(x) -> U(x) + Dr](x) 

Then we can obtain the field equations. Varying the action (0) we have the 



primary equations: 

VM 2 - m 2 (iYD fl + m)ip = (V 2 + m 2 )ip (11) 

VM 2 - m 2 (i7 M V M - m)tp = (V 2 + m 2 )ip (12) 

n(A,-d^)=fJ (13) 

n&>{A lt -d lt <!>) = nu-&>fi (14) 

n<s> = u (15) 
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where we have used the following notation: 



fJ = (16) 



# = - VM2 _ m2 (yV„y + - vyp4 - KM ~ 

- eh^ + erw^ (17) 

Here w^ v is an arbitrary antisymmetric tensor and the last term in the right- 
hand side of eq. ([17]) expresses the arbitrariness of the dependence of the 
current j* from its divergence (note that only d^j® enter the field equations 
- eq- (0)). 

As we mentioned above there are some auxiliary fields in our model. It 
is necessary to exclude them and obtain the final form of the field equations. 
After simple calculations we have from eqs. (O) and (O) : 



(V 2 + m 2 )(z7^V M + m)(p + (M 2 - m 2 )(i'fD^ + m)(p = (18) 
and from eqs. fll4|), ( |15| ) - 



nd^A^ - d^) = n 2 $ - d*j* (19) 

Furthermore from eqs. (p~3"D, ( |I9"D we can obtain 

□ 2 * = 0"(#+#) (20) 
It is easy to verify that the right hand side term in the last equation is zero, 

= (21) 

as a result from field equations fllTD and (|l"2"D. The calculations leading to 
this result one can see in the Appendix. Then 

□ 2 $ = (22) 

and instead of eq. (|T9| ) we have 

Ud^A,-d^) = -d^l (23) 
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Now we can point out the final equations. These are the equations (|IBD, (|T3|) 
and (|22|) for the basic fields and one of the eqs. (|ll]) or (12) for the auxiliary 



field i/j (we ignore the field U). Certainly we must consider our current j 



with excluded ip. The rest of the equations ( pi]) and fl2"3|) follow from the 
basic ones. 

To complete the description of our model, let us introduce new basic fields 
using the above mentioned partial gauge invariance. This we make with the 
help of the following expressions: 

<p(x) = expie$(x).0(x); A^x) = A^x) + <9 M $ (24) 



Then the operator V M turns into in all equations and instead of (|T^ 
and (^3h we have: 



(□ + M 2 )(i7 M 9 A1 + m)4> + erf" Apt = (25) 



M 2 -m 2 



D A = fi (26) 
U&Ap = d% 1 (27) 
respectively. As we can see we have excluded the current (|T7j) from eq. (53) 



using the identity fl2lD . Then equation fl27|) becomes a direct consequence of 



3 Quantization of the model 

We are going to quantize the model described above with the help of the 
perturbation theory in Dirac (interaction) representation f\ For this matter 
we will apply the approach proposed by Bogolubov ||. As is well known, this 
approach involves the formulation of the quantization of the corresponding 
free field theory. In our case we have from eq.(|3]) the following free spinor 
Lagrangean: 

L°(0) = - = =L ^ (d^d^ + d^<j>) + + m)0 + 

2 



y/M' - m 



— 777 — — 

+ — = {M + W>) (28) 

V M 2 — m l 



'Sometimes this representation is called Tomonaga-Shwinger representation. 
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Remark. Let us remind that when we are passing to the quantum theory 
all products of the fields in the Lalgrangeans must be understood as normal 
ones. This convention makes unnecessary the use of any additional symbols 
for the normal product. 

The Lagrangean for the field is the well known Lagrangean for the 
massless vector field and its quantization as electromagnetic field in the 
Lorentz gauge is given, e.g., in p. The field $ satisfying the eq. (^) 
played an important role in the conformal invariant QED and its quantum 
theory is well known too (see for example the works ||). Moreover, we have 
seen that this field can be excluded from our model with the help of the 
gauge transformations (p4|) . On the other hand let us remind that the field 
ijj(x) is an auxiliary spinor field which has to be excluded too. However, the 
Lagrangean (p8[) is more suitable for canonical quantization (containing first 
derivatives only) than one from (]£]). That is why here we will consider the 
theory with the Lagrangean ( p8|) only. The free field equations in our case 
have the form: 

1 (i^dp + m)4> = (□ + m 2 )ip (29) 



m 



2 



2 _ m 2 



{yfdn - m)ip = (□ + m 2 )0 (30) 



The corresponding canonical momenta are 

**(a0 = ^/ M 2_ m 2 d ^ + = M X ))V ( 31 ) 

= VM 2_ m 2 d ^ + ^» = W X ))*-f ( 32 ) 

The anticommutator of the field <f)(x) we denote as 

T a p{x) = {<p a (x),M0)} (33) 

({...} means anticommutator). From equations (^9|) and (|30|) we have the 
following equation for the function r a p(x): 

(□ + M 2 ) (iYd„ + m)T(x) = (34) 
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Furthermore using the canonical commutation relations: 

W(»> <P{y)}x =y Q = {ni>(x),1p(y)}x =y = iS 3 (x)5 al3 (35) 

we have the following initial conditions for the function T 

r<xp(x) U=o = d T al3 (x)\ Xo= o = 
d 2 T al3 (x)\ Xo=0 = (M 2 - m 2 ) 7 o^ 3 (x) (36) 

To obtain the latter we have used the anticommutation relations between the 
fields ip and 0, and the equations they satisfy. These calculations are very 
simple and we omit them here. 

Now it is easy to obtain our function F: 

T aP (x) = -t{tYd„ - m) aP [D m (x) - D M {x)\ (37) 

where D m and Dm are the scalar Pauli- Jordan functions with masses m and 
M respectively. Analogously one can obtain the rest of the singular functions 
from which we will write down here the causal Green function only: 

V c ap {x) = -z(z 7 ^ - m) aP [D c m {x) - D M {x)\ (38) 

D c m and D C M are corresponding scalar causal Green functions for the masses 
m and M: 

Our causal Green function (^) is normalised in such a way that it coincides 
with the propagation function: 

<T0 Q (x),0 /3 (O)> o =r c a/3 (x) (39) 

The formulae (^) and ([38]) show us that the free quantum theory of our field 
is containing indefinite metric in state space. This fact becomes more clear 
if we consider the concrete form of the general solution for the field . It 
has the form 
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<j> a {x) = + 



+ /^-'(PK-iap^t 



+ |/^'(pK + '(«,p) + 



•pi 



ft = Vp 2 + M 2 

p = (pi,P2,p 3 ) P 2 = \jv\ + P 2 2 +: 

Here <t>° a (x) is the usual free spinor field satisfying Dirac equation 

(i^dn + m)(f) (x) = 
and the following anticommutation relation 

{<f>° a (x), $(y)} = -i(ird, - mUD m (x - y) (41) 

The mode-vectors t> ±Q (f2, p) fulfill the usual spinor identities which in our 
case have the form 

(/ P|tT M)^(n,p)=0 (42) 

Moreover this mode- vectors satisfy the normalization condition which can be 
written down in the following two forms: 

E(^ r (n,p))*^'(n,p) = 5„ (43) 

a=l 

or 

4 M 

Y,(vt r (n,p))v± s (n,p) = ±-6 rs (44) 
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The relation 

E(«f-(n,±p))*^(n, T p) = o (45) 

a=l 

is a condition for the orthogonality of the mode-vectors and finally the rela- 
tion 

E « ±r (n,p)(t^-(n,p)) = 7P ; o (46) 

r=l,2 Zii 

expresses the completeness of these vectors. 

The quantities a +r (p), b +r (p), c~ r (p) and (i~ r (p) entering the general 



solution ( pEOP are operators in the quantum theory and they must satisfy the 



following anticommutation relations: 

{a +r (p),a- s (q)} = -Q5 rs 5\p - q) 

{6 +r ( P ),6-(q)} = -ri5 rs 5 3 ( P -q) (47) 

{c +r ( P ),c- s (q)} = -^ rs 5 3 (p-q) 

<T s (q)} = -Q5 rs 5 3 (p-q) 

(all other anticommutators are zero) according to the (|33|). In the above 
written formulas we have made the following notation 

a' r (p) = (a+ r (p)Y 

b- r (p) = (b+ r (p)Y (48) 
c +r ( P ) = (c- r (p)Y 

d^(p) = (d- r (p)y 

The new quantities (a +r (p))*, (6 +r (p))*, (c" r (p))*, (d _r (p))* appear in the 
field <f> a (x) Dirac conjugated to the field ( fTO|) and the sign * has the meaning 
of Hermitean conjugation. 

Now let us describe the structure of our state space S. First of all S is 
containing the subspace H which coincides with the Fock space of the free 
quantum field 4>°(x). This is the usual state space of the free quantum spinor 
field with positive metric. Its vacuum state | > is common to the whole 
space S and satisfies the following additional conditions: 



a - r (p) | > = r r (p) | >= c" r (p) | >= cT r (p) I >= (49) 
< | a +r (p) = < | b +r {p) =< | c +r (p) =< | d +r {p) = 
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Then we can see that S contains an additional subspace G in which the basis 
is formed from all monomials of creation operators (a +r (p), 5 +r (p), c +r (p), 
d +r (p)) acting on the vacuum. The corresponding operators with the sign "- 
" are the annihilation operators. According to the our commutation relations 
given above, we can do the following conclusions: 

i) The metric of the space G given with the usual scalar product of the 
Fock space is indefinite. 

ii) The spaces H and G are orthogonal to each other and 

S = H® G 

It is obvious that the vectors of the space S can not be used as physical states 
because of the presence of the ghosts from subspace G. 

However, we can define the physical space to coincide with the subspace 
H and consider only its vectors as edge physical states. Then it is easy to 
verify that 

< ph' | 0(x)7 M 0(x) I ph >=< pti | 0°(x)7 / "0°(x) | ph > 

according to the eq. (0) and of course 

df, < ph' | <p(x)Y<P(x) \ph>=0 (50) 

where | ph > and | ph' > are two arbitrary states from the physical space 
H. Moreover, in this case the relation (pT|) define the space H and can be 
considered as a defining condition for the physical states. 

Passing to our model with interaction it is naturally to generalize the 
condition ( |50"D for the definition of the physical states in the interaction case. 
This we can do defining the latter with the help of the relation, formally 
coinciding with relation floUD but with <p(x) satisfying the equations 
([2~?D . Such definition is in accordance, except for the free case, with the 
quantum theory in the interaction representation, where the spinor physical 
space coincides with the space H. 

As we mentioned above we consider our model only on the physical space. 
This means that all matrix elements of operators and their products having 
some physical meaning must by taken between physical states only. Then 
the nonphysical states such as ghost ones will give contribution to the inter- 
mediate virtual states only. 
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Now we are going to see what happens with our model in the physical 



space. According to the condition (50) we can obtain there the form of the 
equations Q2T5D and fl2"T|): 



< pti | UA^x) - j*(x) \ph>=0 (51) 

and 

< pti | nd^A^x) \ph>=0 (52) 

The last equation shows us that in the chosen physical space automatically 
appears the gauge condition such as 

na^Am = 

This condition contains the Lorentz gauge fixing 

d^Afj, = 

which means that the our physical space has a subspace in which equation 
(|5"ID coincides with the Maxwell equation in the Lorentz gauge. The vectors 
belonging to this restricted physical space are denoted here as | pho >. Then 
instead of eqs. (|5"ID and (|5^) we have 

< pK I n A( x ) -tf( x ) I P h o >= 

=< pti | d v F v „{x) - fj(x) | ph >= (53) 
F = d A — d A 

and 

<9 M < pti | Ap(x) | ph >= (54) 

The construction described here obtains its concrete form in the interac- 
tion representation, where, as is well known, the Lagrangean and the other 
physical quantities are expressed through the corresponding free quantum 
fields in the Heisenberg representation. In this case the mentioned above 



restricted physical space Ho in which equations fl5H|) and (|54| ) take place, 
coincides with the following space with nonnegative metric: 

H = H®R a (55) 
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where the space R a can be defined as 

N 

Ra=®®l\ (56) 
jV>0 k=0 

iS =| it >=| o > 

and where lj. are the one particle spaces for the field A^x). The arbitrary 
one particle state in this case has the form 

| I 1 >= J A+{x)F{x)d 4 x | > (57) 

where A+(x) is the positive frequency part of the field A^(x) and l^(x) are 
arbitrary test functions fulfilling the condition: 

d^F(x) = (58) 

i.e., 

d»A-(x) | V >= - / d»D~{x - y)Uy) d ^ = 

% J 

Here A~(x) and Dq(x — i/) are the negative frequency parts of the field Ap(x) 
and the massless Pauli Jordan function respectively. 

Our physical space is the subspace of the whole quantum state space Q 
which in the considered representation has the structure 

Q = S®F A ; H CQ (59) 

where Fa is the Fock space of the quantum vector field A^(x) satisfying the 
d'Alambert equation. 

Now we can formulate the main result of the present paper: 
The quantum model with the action (^|) and Lagrangeans (|3|), (|J) and 
fllQl) , leading to the field equations with the higher derivatives ( p5|) , (|26|) and 
(P?D, contains in its own state space Q the subspace H with nonnegative 
metric in which our model coincides with the spinor QED. From this point 
of view the latter is immersed in our model described above. It is necessary to 
mention that this QED is in regularized form, because the spinor propagation 
function is given by expression (|38| ) and, as is easy to see, it coincides with 
the Pauli- Villars regularized one. That is why the Feynman diagrams are 
free from divergences in our case. In the next section we are going to discuss 
these questions. 
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4 Regularized QED 



In this section we would like to compare the described above regularized 
QED (RQED) which is immersed in our model, with the standard one. This 
comparison we are going to do on the level of Feynman diagrams and it helps 
to understand better our reasons to consider such theory. 

As we have seen, the physical space of RQED in the interaction represen- 
tation is composed from Fock spaces of the free spinor <fi° a (x) and free photon 
A^(x) fields in Lorentz gauge. This space is with nonnegative metric and 
gives us the possibility to do the next steps to obtain the physical space of 
the free electrons, positrons and photons, e.g., passage to the Coulomb gauge 
fixing, factorization of the zero norm state vectors and so on (see ref. 0); 
the steps usually done in standard QED too. Then it is easy to notice, that 
our physical space is the same as in standard QED. There is no difference 
also between the photon propagation functions in the two theories. In our 
case this function is 

= V^D c (x) 

where rj^ is the Minkowski metric tensor and D^(x) is the causal Green 
function of the d'Alambert equation. 

Remark. In our model we have used the Lorentz gauge fixing as the most 
simple. However, there are no obstacles to use other forms of the photon part 
of our RQED. For instance if we had used, instead of Lagrangean (f|), a new 
one as follows: 

L V (A) = -1/2&*(A V - dvtyd^A" - d v $) + ~0"(A„ - d u <b)d»{A, - dfl) 

the corresponding equations for the photons on the physical states would 
have obtained the form 

with the following gauge fixing condition 

= 

This form of the equation for the electromagnetic field A^ is well known too 
(see ref. [ffDI) 
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The main difference between RQED and the standard QED is in the form 
of the spinor propagation function. As it is seen from eq. (|38D the form of 
our spinor propagation function coincides with the well known Dirac one but 
regularized with the help of the Pauli-Villars regularization without perform- 
ing the limit M — > oo. In the momentum representation the function (|38|) 
has the behaviour p~ 3 for p — > oo, so RQED is free of ultraviolet divergences. 
However, a new parameter M appears here. For the understanding of the in- 
fluence of M over the RQED we are going to consider all Feynman diagrams 
arranged in two groups. In the first group we put all diagrams which remain 
finite in the limit M — > oo. These diagrams correspond to converging ones in 
the standard QED. Here it is possible to choose M large enough so that the 
corresponding diagrams from the two QED's will coincide with each other 
on such a range of the values of the external momenta, which could have 
any experimental significance. That is why we can say that the first group 
of diagrams RQED differs from standard QED only for very big external 
momenta, the range of values of which is defined by the value of the mass 
parameter M. 

In the second group of Feynman diagrams we put those which correspond 
to the divergent diagrams of QED. These diagrams increase unlimitedly for 
large M. In the standard QED over these diagrams one applies the renor- 
malization procedure. In the our case, because of the finiteness there is no 
similar necessity in the RQED. That is why all nonuniqueness appearing 
in QED after the infinite renormalization is described here through the pa- 
rameter M only. In our opinion this is the main difference between the two 
considered here QED's. We would like to demonstrate this difference with the 
help of some example. For this we choose the second order electron-positron 
self-energy Feynman diagram, i.e., one of the basic divergent diagrams in 
standard QED. 

As usual the mentioned diagram we denote as E aj g(p), where p is the 
external electron-positron momentum. Without giving here the calculation 
which is well known, we will start from the following expression for this 
quantity in our case: 



The integration in the right hand side of the last equation can be taken and 
for us is interesting the result when M 2 3> p 2 , i.e., in the low energy range. 




(60) 
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Then we have 

e 2 M 2 

= 77—2^ m - 7%)^— + tf(p) (61) 
1671^ rrr 

where we have denoted by the part which does not increase with M 

£'(p)= / de(2m-^)fn 2 2 +Q(tt) (62) 
87T 2 jo m z — £;r m 

The corresponding diagram in the standard QED after renormalization has 
the form 

E(p) = Cl (Y*p„ - m) + c 2 + S'(p) (63) 

where S'(p) is the same as in eq. (|62|) (see ref. 0). Here C\ and c 2 are 
arbitrary finite constants in the result of the renormalization. Comparing 
eq. ( |63|) with eq. ([H]) we can see that in our case both constants depend on 
M as follows 

-e 2 , M 2 3e 2 m, M 2 

C\ = Tiln — - c 2 = Trln — - 

I67H m 1 I67H m 2 

This means that in RQED we have no nonuniqueness appearing in the process 
of the calculation of such a diagram. 

Remark: It is possible to do finite renormalization of our model. Then 
the analogous nonuniqueness will appear in RQED too. However, the finite- 
ness of the considered diagrams makes such procedure unnecessary. 

The only quantities which can be considered as arbitrary parameters are 
the mass parameters M and m. Recalling that E Qi g(p) is the second order 
correction to the free spinor operator we can write down the corresponding 
corrected Green function in the momentum representation: 

G{P) = j^{v 2 -M 2 ){YV»-m)-nv) (64) 

Now we can calculate the pole-point of this function which will represent the 
electron-positron physical mass m p /j. Up to second order, this mass can be 
obtained in the form: 

m ph = m + $ m 

3e 2 m 7 M 2 5e 2 m 

6m = — 77TT — 2 ~ ~ooT 65 
1671^ m z S2tt z 
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It is well known that the point p 2 = m 2 is a branching point of 
and because of the absence of the renormalization it is impossible to set 
rriph = m. That is why there are bare m and dressed m p h electron-positron 
masses in RQED, connected by expression (|65|). The dressed mass m ph is a 
pole of the Green function (Q) and it is less than the bare mass m which 
is a branching-point of the same Green function. Here these points are au- 
tomatically different and there is no necessity in RQED to introduce small 
photon mass to reach this difference, as it takes place in the standard QED. 
Perhaps this is the most significant property. If we reformulate it, this prop- 
erty means that there are no infrared divergences in the considered diagram 
(see for this ref.0). 

5 Conclusion 

We defined a model with higher derivatives. Our analysis shows that in the 
suitable chosen physical subspace of states, this model coincides with QED. 
The only difference between QED and our model (calling here RQED) is the 
usage of the regularized causal spinor Green function in the latter. However, 
in our theory we can not take M — > oo, because the mass parameter M enter 
the initial spinor part of the Lagrangean and it has a finite value. After simple 
analysis we saw that RQED is free from ultraviolet divergences. Choosing 
the parameter M sufficiently large, we can say that RQED and QED have 
the same behaviour in the low energy range. We are hopeful that RQED will 
have a behaviour in the high energy range better than QED, because of the 
faster decrease of our spinor propagator for large momentum. 

On the example of the spinor selfenergy second order Feynman diagram 
we saw that there is no infrared divergence too. 

From the obtained results we can see that there exists an alternative 
way to consider given particle theory in which the field operator describes 
not only the physical one. Then the physical theory turns out as immersed 
into wider theory, such that the latter has contribution to the intermediate 
(virtual) states only. We have demonstrated here this approach on the well 
known theory for comparison. However, there are theories such as nonrenor- 
malizable ones, where this approach could turn out to be the only possible. 
For example, applying our propagator to the four-fermions interaction, the 
corresponding theory will have no divergences at all. 
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APPENDIX 

First of all let us write down the full current jff + j * in the terms of fields 
defined in eq. (0). Then we have 

it + ft = -ed>>h»<t> ~ e&bvZ- 
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VM 2 - m 2 
where by £(x) we have denoted 

£(x) = exp[—ieQ(x)].ip(x) 

Then the equations (pi]) and (0) can be written down as follows: 

VM 2 - m 2 (iYd IM + m)(f> = (□ + m 2 )£ - eA^cj) (Al) 



VM 2 - m 2 {i-f fl d fl - m)£ = (□ + m 2 )cf> (A2) 



We can do the verification of eq. (|21|) with direct calculations. For this let 
us write down the Dirac conjugates to eqs. (A.l) and (A. 2) 



VM 2 - m 2 {-id ll <j)Y + m<p) = (□ + m 2 )£ - e^A^ (A3) 



VM 2 -m 2 (-id^Y - m£) = (□ + m 2 )0 (A4) 

Now we must multiply eqs. (A.l) and (A. 3) by from left and — from 
right respectively and then add the results. Analogously we must proceed 
with eqs. (A. 2) and (A. 4). It is easy to see then, that the sum of the two 
obtained in this manner identities coincides with the expression: 
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